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Abstract 

Proof of factorization of soft and collinear divergences in non-equilibrium QCD may be necessary 
to study hadronic signatures of quark-gluon plasma at RHIC and LHC. In this paper we prove 
factorization of soft and collinear divergences in non-equilibrium QED by using Schwinger-Keldysh 
closed-time path integral formalism in the background field method in pure gauge. 
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I. INTRODUCTION 



RHIC and LHC heavy-ion colliders are the best facilities to study quark-gluon plasma 
in the laboratory. Since two nuclei travel almost at speed of light, the QCD matter formed 
at RHIC and LHC may be in non-equilibrium. In order to make meaningful comparison 
of the theory with the experimental data on hadron production, it may be necessary to 
study nonequilibrium-nonperturbative QCD at RHIC and LHC. This, however, is a difficult 
problem. 

Non-equilibrium quantum field theory can be studied by using Schwinger-Keldysh closed- 
time path (CTP) formalism However, implementing CTP in non-equilibrium at 
RHIC and LHC is a very difficult problem, especially due to the presence of gluons in 
non-equilibrium and hadronization etc. Recently, one-loop resummed gluon propagator in 
non-equilibrium in covariant gauge is derived in |3|, |4 ] . 

High pt hadron production at high energy e + e~, ep and pp colliders is studied by using 
Collins-Soper fragmentation function [5, (], J For a high px parton fragmenting to hadron, 
Collins-Soper derived an expression for the fragmentation function based on field theory and 
factorization properties in QCD at high energy. This fragmentation function is universal 
in the sense that, once its value is determined from one experiment it explains the data at 
other experiments js]. 

Recently we have derived parton-to-hadron fragmentation function in non-equilibrium 
QCD by using Schwinger-Keldysh closed-time path integral formalism \$. This can be 
relevant at RHIC and LHC heavy-ion colliders to study hadron production from quark- 
gluon plasma. We have considered a high p? parton in QCD medium at initial time r with 
arbitrary non-equilibrium (non-isotropic) distribution function f(p), fragmenting to hadron. 
The special case f(p) = -tw — corresponds to the finite temperature QCD in equilibrium. 

e T ±1 

In order to study hadron production from quark-gluon plasma using the non-equilibrium 
fragmentation function we will need to prove factorization of fragmentation function in non- 
equilibrium QCD. Factorization refers to separation of short distance from long distance 
effects in field theory. In this paper we will prove the factorization of soft and collinear 
divergences in non-equilibrium QED. The proof of factorization of fragmentation function 
in non-equilibrium QCD at RHIC and LHC will be addressed elsewhere. 

The paper is organized as follows. In section II we briefly describe factorization of soft 
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and collinear divergences in QED in vacuum. In section III we describe Schwinger-Keldysh 
closed-time path integral formalism in non-equilibrium quantum field theory relevant for 
our purpose. In section IV we prove factorization of soft and collinear divergences in non- 
equilibrium QED by using Schwinger-Keldysh closed-time path formalism in the background 
field method in pure gauge. Section V contains conclusions. 



II. FACTORIZATION OF SOFT AND COLLINEAR DIVERGENCES IN QED IN 
VACUUM 

The Wilson line is given by 0, Q 

$[x»] = V exp[-ie f° dX h ■ A(x^ + WX)] (1) 

J — oo 

where is a x^ independent four vector. Eq. (pQ) can be written as 

= V exp[-ze [ dXh- A{x" + h^X)} 

J — oo 

= V exp[-ie [° dXh- e Xh9 A(x»)] = V exp[-ie-^—h ■ A(x% (2) 
J— oo ri • o 

Since the vector is free we can choose it to correspond to various physical situations. For 

example if we choose = n^, where is a fixed lightlike vector, we find from eq. ([2]) the 

phase factor 



lu(x) = -n ■ A(x). (3) 

n ■ o 



Using the Fourier transformation 



we find from eq. 



V = eu(k)=ie rA^(k). (5) 

Th ' fx. 



Note that eq. ([5]) is precisely the eikonal vertex for a soft photon with momentum k in- 
teracting with a high energy electron jet moving along the direction . In the soft 
photon approximation in [6( the fixed lightlike vector is taken to be having only "+" or "-" 
component: 

n» = (n + ,n-,n T ) = (1,0,0) or = (n + , n~, n T ) = (0, 1, 0). (6) 



Now we will show that when the classical background field A^(x) is a pure gauge given 

by 

A„(x) = d„cu(x) (7) 

we can reproduce eqs. ([3]) and ([5]) which appear in the Wilson line eq. ([2]). Now multiplying 
a four vector /i M from left in eq. ([7]) we find 

h-A(x) = h-duj(x). (8) 

Dividing h ■ d from left in the above equation we find 

uj(x) = -—h-A(x). (9) 
h ■ o 

Since the four vector is free we can choose it such way that the pure gauge can represent 
Feynman rules involving soft photons or collinear photons. For example, when = 
where n M is a fixed vector (see eq. (J6])), we find from the above equation 

uj(x) = ^-gK ■ A(x) (10) 

which reproduces eq. ([3]) which appears inside the Wilson line in eq. ([2]). Hence we have 
established the correspondence between the Wilson line and the classical background field 
A^x) in pure gauge in the context of soft divergences. 

Similarly, if we choose = n 1 ^, where is a non-light like vector 

n>B = (n£,ni,0), (11) 

we reproduce the Feynman rules for the collinear divergences {(J. This establishes the 
correspondence between the Wilson line and the classical background field A^(x) in pure 
gauge in the context of collinear divergences. 

The generating functional in QED in the presence of background field A^x) is given by 

H 

Z[A, J, V,V} = J [dQ] [#][#] e i / d4 " [ -3^[Q]-^^^) 2 +^[A+Q^+J-Q+^+^] (12) 

where 



Ef[A + Q] = (i$-eQ-e4). 
4 



(13) 



For the purpose of studying soft and collinear divergences it is understood that the integra- 
tion in Z[A, J,i],fj\ is only over those functions of the quantum photon field Q(x), electron 
field ip(x) and positron field ip{x) such that their Fourier transforms Q(k), ip{k) and ip(k) 
vanish in the soft region. Under the following transformation of the fermion fields and the 
sources 

we find from eq. (TT2I) 

Z[A,J,rf,ff] = Z[J,r),fj\. (15) 

The correlation function in QED in the presence of background field A^{x) is given by 
6 6 

x -, x x i \ Z [ A , J, V, fj] \j =r] =fi=0 =< ^(X 2 )^(X 1 ) > A . (16) 

Similarly the correlation function in QED (without the background field) is given by 
S S 

5^2) H^) =< >A =° ■ (17) 

Hence we find from eqs. ( JTBl . (TIBl) and (TlTl) 

< V(ar 2 )^(xi) >a = e~ ie ^ < ^(z 2 )^(xi) >a=o e ie ^ (18) 

where all the A^ix) dependence has been factored into ^- %ew ^) and e ieu) ^ xi \ Using eqs. ([9]) 
and (j2D we find 

r° 

< Tp(x 2 )ip(xi) > A = exp[-ze / dX h ■ A(x% + h^X)] 

J —oo 

' I 





x < ip(x 2 )ip(xi) >a=o x exp[-ze / dX h ■ A{x\ + h v X)\. (19) 
This proves factorization of soft and collinear divergences in QED lOj. 



III. SCHWINGER-KELDYSH CLOSED-TIME PATH INTEGRAL FORMALISM 
IN NON-EQUILIBRIUM QUANTUM FIELD THEORY 

Unlike pp collisions, the ground state at RHIC and LHC heavy-ion collisions (due to the 
presence of a QCD medium at initial time t = ti n (say ti n =0) is not a vacuum state |0 > any 
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more. We denote [in > as the initial state of the non-equilibrium QCD medium at ij n . For 
example, if the system at RHIC and LHC at initial time is space translational invariance 
the non-equilibrium distribution function f(k) of a parton (quark or gluon), corresponding 
to such initial state can be written as 

< a\k)a(k') >=< in\a\k)a(k')\in >= f {k) (2n) d - 1 5 {d - 1} (k - k'). (20) 

Finite temperature field theory formulation is a special case of this when f(k) = — . 

eT ±1 

Consider the time evolution of the density matrix 

dp 

= [Hi, Pi p(-oo) = Po (21) 

where Hi is the interaction hamiltonian. The formal solution is 

p(t) = S(t, -oo)p S(-oo,t), S(t, -oo) = T exp[-i f dt' H^t')}. (22) 

J — oo 

The density matrix p is in interaction picture. The average value of an operator L in the 
interaction picture is given by 

< L(t) >= Tr[p{t)L{t)i i^- = [H , L(t)] (23) 

where Hq is the free hamiltonian. Since in many situations we deal with correlation function 
of several fields at different times, it is useful to transfer all the time dependence to field 
operators and consider the density operator as independent of time, i. e. to go to Heisenberg 
representation. For the time independence of the density matrix we can take the value 
of the matrix determined by expression eq. (|2"2"|) at a certain fixed instant of time, for 
example, t = t in = 0, having thus included in it all the changes which the distribution p 
had undergone when the external field and the interaction in the system were switched on. 

Consider the medium average of T— products of several operators. Using the Heisenberg 
density matrix one finds 

< T[L(t)M(t')...] >= Tr[pT[L(t)M(t') ....]] = Tr[S(0, -oo)p S(-oo, 0)T[L(t)M(t' )....}} 
= Tr[p S(-oo, 0)T[L(t)M(t')....}S(0, -oo)]. (24) 

Going over to the operators in the interaction picture 

< T[L(t)M(t')...} >= Tr[p S(-oo, 0)T[S(0, t)Lj(t)S(t, t r )Mi(t')....]S(0, -oo)} 

= Tr[p T c [5 c L / (t)M / (t / )--]] (25) 



where T c is the complete contour from 

- oo -> t -> t'.... -> -oo (26) 

and S'c is the complete 5— matrix defined along T c . 

To deal with Feynman diagram and Wick theorem it is useful to split the time interval 
to "+" and "-" contour; where "+" time branch is from — oo to +00 where (time) T— order 
product apply and "-" time branch is from +00 to —00 where (anti-time)T— order product 
apply. 



A. Generating Functional in Non-Equilibrium Scalar Field Theory 

Consider scalar field theory first. Since there are two time branches there are two fields 
and two sources and hence four Green's functions. Let us denote the field <f>+(x) and the 
source J+(x) in the "+" time branch and <f>-(x) and J-(x) in the "-" time branch. The 
generating functional is given by 

Z[J+, J-,p] = J \d<p + }[d<p_} exp[i[S[(p + ] - S[<f>-] + J d 4 xJ + <p + - J d 4 xJ„<p„}] < 0+,O|p|O,0_ > 

(27) 

where S[<f>] is the full action in scaler field theory and \(f>±,0 > is the quantum state corre- 
sponding to the field configuration <fi±(x,t = 0). 

In the CTP formalism in non-equilibrium there are four Green's functions 

G ++ {x,x') = .ff^ + \ J " , . P] / . =< in\T<t>{x)4>{i!)\in >=< T<p(x)^x') > 
i z oJ + {x)J + {x') 

G-(x,x') = S 4ttTyA n =< MTm<t>^')\in >=< T4>{x)4>{x>) > 
{—t) z dJ-{x)J_{x') 

G + -(x,x') = tf/Vw P ( a =< ™\<KJ)<K*)\in >=< <t>{x')4>{x) > 
—t z dJ + {x)J-{x l ) 

G -^ x ') = S SrrVTT a =< M<t>{x)4>{x')\in >=< <t>{x)<t>{x') > (28) 
— i A oJ_{x) J+yx 1 ) 

where T is the time order product and T is the anti-time order product given by 

T(f)(x)(f)(x') = 6{t - t')(f)(x)<f)(x') + 6{t' - t)(f)(x')(f)(x) 

T<t){x)<t){x') = 9{t' - t)(f)(x)(f)(x') + 6{t - t')<f){x')<f){x). (29) 
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B. Generating Functional in Non-Equilibrium QED 



Extending the same analysis to QED we find the generating functional in non-equilibrium 
QED 

Z[J+,J_, 77 + ,77_,77 + ,?7_] =/'[dQ + ][dQ_][d^ + ][#_][d'0+][#-] x 

exp[* / d*x[-~(F^[Q + ] - F%[Q_\) - ^(( W) a - (3 M Q-") 2 ) + 
-$_i?[Q_]V>- + J+-Q+- fj + ^+ - + r] + ^ + - 77_^0_]] 

x < Q + , V+, 0| p |0, > ■ (30) 

where Q M is the photon field and J M is the corresponding source, t/> is the electron field and 77 
is the corresponding source, p is the initial density of state. The state \Q±, i()±, V>±, > is the 
quantum state corresponding to the field configurations Q^(x, t = t in = 0), -^(x, t = t in = 0) 
and ip(x,t = ti n = 0) respectively. 

IV. FACTORIZATION OF SOFT AND COLLINEAR DIVERGENCES IN NON- 
EQUILIBRIUM QED 

The generating functional in the background field method of QED is given by eq. ( fT2l) . 
Extending this to non-equilibrium QED we find 

Z[p,A,J + ,J_,i] + ,r]_,f] + ,f]_] = j[dQ + )[dQj\[d$ + )[dTp^)[di) + )[d'ip„) x 

exp[* J d i x[~{F%[Q + ] ~ F%[Q-]) ~ ^(( W) 2 - (d,Q-n 2 ) + + 
-$-If[Q- + A-]tl>- + J+-Q+- J--Q- + v+ij + - fj-il>- + V+4>+ - V-*P-}} 
x < Q + , if> + , tp + , 0| p |0, Q- > (31) 

where p is the initial density of states. The fermion fields and the sources transform as 

= U- ip-, = UZ\ rj'_ = U- ff_ = fj_ UZ 1 (32) 

where 

U + = e - ieuJ+{x \ U- = e~ ieaJ - (x) . (33) 
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Changing ijj± — > ijj' ±} ip± — > ip' ± but keeping Q± fixed we find from eq. (l3Tj) 
Z[p,A, J+, J-,r] + ,r}^,fj + ,fj_) = J[dQ + )[dQ^) [#+][#-] [#+][#_] x 

exp[* / d 4 x[-i(F^[g + ] - F^QA) - ^((^Q+^) 2 - (9,Q^) 2 ) + ^ + Z/[Q + ]^ + - $-Lf[Q-\1>- 
+J+ ■ Q + - J_ ■ Q_ + 77+^ - 77_V- + 77+^ - 77^]] 

x < Q +7 iP' + ,fi + ,0\ p \0,fi_,iP'_,Q_ > . (34) 
Note that the state |0, Q_ > is at initial time t = t,i n = 0. We find 

Z[p, A,J + ,J-,rj + ,r]-,f} + ,f}-} = J x 

exp[* J d*x[-±(F%[Q + ] - F%[Q.]) ~ ~ (^Q- M ) 2 ) + $+I?[Q+]il)+ ~ $-Lf[Q-\$- 

+J+ ■ Q + - J_ ■ Q_ + f] + f + - fj-.iI/_ + ~ V-ft-]] 

x < Q + ,^ + ,0| p |0,V>_,^_,g_ > (35) 
which gives by using eq. fl32l) 

= |[dg+][dg_] [#+][#_] [#+][#_] x 

ex P [* J d*x[-±(F%[Q+] - FlAQA) - ^((<V?+ M ) 2 - (d,Q-n 2 ) + ^[Q + ]^ + ~^[Q~^- 
+J+ -Q+- J--Q-+ fj+ip+ - f)-ip- + - v-$-\] 

x < Q + ,^ + ,V + ,0| p |0,^_,V-,Q- > (36) 

Hence we find from eq. f l36l) 

Z[p, A, J + , J_, 77+, 77^, 77+, 77^] = Z[p, J + , J-,r} + , ?7_, 77+, (37) 

in non-equilibrium QED in the presence of background field in pure gauge. 

The correlation functions in non-equilibrium QED in the presence of background field 
Apia) is given by 
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z[p,AJ+,J-,v+,v-,v+,v- 



J + =J-=r) + =T)_=r) + =ri-=0 



5r] r (x 2 ) drjstx! 
= < i> r (x 2 )i>s(xi) >A = < in\ip r (x 2 )ip s (x 1 )\in> A 

(38) 

where r,s = +, — are the closed-time path indices. The correlation functions in non- 
equilibrium QED is given by 
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Z[p,J + ,J_,r} + ,r}„,r] + ,r}„ 



\J+=J- =T7 + =»?_ =-q + =?7_ =0 



Sr) r (x 2 ) Sl] 3 (Xi) 

= < ip T (x 2 )i' s (x 1 ) > A =o = < in\ip r (x 2 )'ips(xi)\in >a=o ■ (39) 



Hence we find from eqs. (|37j) . (!38|) and (|39|) 

<m|^r(x 2 )^(xi)|m> A = e - ie ^ (x2) <m|^ r (x 2 )^ s (x 1 )|m>^ =0 e ieLJfl(a;i) . (40) 

where all the A^(x) dependence has been factored into e~ teuJr( - X2 ^ and e ieu)s<yXl \ Using eqs. 
© and © we find 

< in\ip r (x2)ips(xi)\in >a = exp[— ie / dX h ■ A r (x2 + h^X)} x 

J — oo 

r — oo 

< in\vp r (x2)ip s (xi)\in >a=o x exp[— ie / dX h ■ A s (x\ + h u X)} (41) 

Jo 

which proves factorization of soft and collinear divergences in non-equilibrium QED. This 
is similar to eq. ([1^]) except that the closed-time path indices r, s = +, — have appeared in 
eq. (|4"T|) . Note that the repeated indices r and s in eq. (|4"T|) are not summed. 



V. CONCLUSIONS 

In order to study hadron production from quark-gluon plasma using the non-equilibrium 
fragmentation function [9] we will need to prove factorization in non-equilibrium QCD. 
Factorization refers to separation of short distance from long distance effects in field theory. 
However, before proving factorization in non-equilibrium QCD we need to show a similar 
factorization in non-equilibrium QED. In this paper we have proved factorization of soft and 
collinear divergences in non-equilibrium QED by using Schwinger-Keldysh closed-time path 
integral formalism in the background field method of QED in pure gauge. 



The proof of factorization of fragmentation function in non-equilibrium QCD at RH 



and LHC will be addressed e 



from quark-gluon plasma llj at RHIC and LHC. 



sewhere. This may be relevant to study hadron production 
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